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MULTI-CHANNEL FREQUENCY DOMAIN EQUALIZER 
FOR RADIO FREQUENCY INTERFERENCE CANCELLATION 



Background of the Invention 

1. Field of the Invention 

This invention relates generally to radio frequency interference impairments 
associated with digital subscriber lines, and more particularly to a method for exploiting 
deterministic spreading of the fast Fourier transform (FFT) at the receiver to cancel 
correlated noise sources such as RFI. 

2. Description of the Prior Art 

ITU Telecommunications Standardization Sector, G.tesLbis., Study Group 4, 
Question 15, 1999, has proposed various radio frequency interference (RFI) impairments 
for incorporation into digital subscriber line (DSL) performance tests. For asymmetric 
DSL (ADSL) systems such as described in ITU Telecommunications Standardization 
Sector, G.992J: Asymmetrical Digital Subscriber Line (ADSL) Transceivers., Study 
Group 4, Question 15, 2000, RFI is assumed to be caused by narrowband amplitude 
modulated (AM) radio stations, which overlap the downstream ADSL band from 540- 
1 104 kHz. Additional RFI sources include, but are not limited to, aliases from AM 
stations operating above 1 104 kHz, amateur (HAM) transmitters, televisions, computer 
equipment, and various other (relatively) low frequency signals (e.g., aeronautical 
radiolocation signals). 

At the ADSL transmitter, data is assembled into subchannels, transformed by an 
inverse fast Fourier transform (IFFT), and a cyclic prefix is added. After propagating 
through a communication channel, the opposite operations are performed by the receiver. 
The cyclic prefix is removed, and the received data is transformed by the fast Fourier 
transform (FFT). If the communication channel is shorter than the cyclic prefix and the 
noise is uncorrected on the individual subchannels, then the transmitted data can be 
recovered by a one-tap equalizer (the frequency-domain equalizer or FEQ), such as 
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described by T. Starr, J. Cioffi and P. Silverman, Understanding Digital Subscriber Line 
Technology, Upper Saddle River, NJ: Prentice Hall, 1999. 

The noise that results from RFI at the receiver is deterministically spread by the 
FFT to neighboring subchannels. As a result, instead of one (or a few) subchannels being 
effected by the narrowband impairment, many subchannels are effected. The resulting 
loss of signal-to-noise ratio (SNR) can unacceptably limit the capacity of the 
communication channel. 

In view of the foregoing, a need exists for a scheme to exploit the deterministic 
spreading of the FFT at the receiver to cancel correlated noise sources such as RFI in 
order to increase the effective data communication rate between subscribers using shared 
media. 
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Summary of the Invention 



The present invention is directed to a method of estimating the multichannel FEQ 
in block and iterative fashions. Specifically, algorithms are developed to estimate the 
multichannel FEQ. The estimated multichannel FEQ provides improvements in SNR 
over that achievable using standard single channel FEQ. 

According to one embodiment, a multichannel FEQ is formulated when there is 
correlation in the noise matrix (e.g., from a correlated noise source such as RFI). 
Equalizers which have nonzero elements in addition to the main diagonal are considered. 
Each row of the multichannel FEQ is recovered via choosing the nonzero columns in 
addition to the column directly correlated with the input through the channel. Columns 
corresponding to subchannels having strong correlations with the column directly 
correlated with the input through the channel are most preferably chosen (generally 
neighboring subchannels, or subchannels corresponding to strong RFI sources). As the 
RFI and signal power become closer, not transmitting data on a subchannel allows that 
subchannel to be used more effectively with other subchannels in the multichannel FEQ. 

As used herein, FEQ means a complex multiply to adjust for the gain and phase of 
the effective channel response in each sub-symbol. The FEQ does not enhance the SNR; 
it only scales and rotates each sub-symbol so that a common decoder can be used. 

As used herein the following notation will have the following meanings. Italics 
will be used for scalars, and boldface for vectors and matrices. A superscript _1 will refer 
to the pseudoinverse, * to conjugation, r to transpose, and H to conjugate transpose. The 
operator diag() forms a diagonal matrix with the arguments being the elements of the 
main diagonal, wherein E[] is the expected value of its argument, and Amax() is the 
maximum eigenvalue of its argument. The term a:b will be used to denote a contiguous 
subset of a vector or dimension of a matrix, and a ":" by itself refers to all of the elements 
along that dimension. For example, if 



c = [c(l)...c(L)] r , 



then, 
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c(:) = c, and 

c(a:b) = [c(a)...c(b)] r 

To select a (possibly) noncontiguous subset of a vector or dimension of a matrix, a vector 
of the desired indices will be used as the argument of that dimension. For example, if 

k = [k(l)...k(M)] T 
is a vector of integers, then 

c(k) = [cWl))...c(t(M))] r 
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Brief Description of the Drawings 

Other aspects and features of the present invention, and many of the attendant 
advantages of the present invention, will be readily appreciated as the same become 
better understood by reference to the following detailed description when considered in 
connection with the accompanying drawing wherein: 

Figure 1 is a graph illustrating a simulated subchannel SNR for additive white 
Gaussian noise (AWGN) with no RFI, AWGN with RFI using single channel FEQ, and 
AWGN with RFI using multichannel FEQ in accordance with one embodiment of the 
present invention; 

Figure 2 is a system diagram depicting a multichannel FEQ integrated with an 
ADSL communication system according to one embodiment of the present invention; and 

Figure 3 is a flow diagram illustrating a method of estimating multichannel FEQ 
coefficients for the n th subchannel, g(«), according to one embodiment of the present 
invention. 

While the above-identified drawing figures exemplify characteristics associated 
with particular embodiments, other embodiments of the present invention are also 
contemplated, as noted in the discussion. In all cases, this disclosure presents illustrated 
embodiments of the present invention by way of representation and not limitation. 
Numerous other modifications and embodiments can be devised by those skilled in the 
art which fall within the scope and spirit of the principles of this invention. 
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Detailed Description of the Preferred Embodiments 

The present invention is best understood by considering a system with n = 1, . . 
N subchannels which transmits data at time t = 1 , . . ., T. If subchannels are arranged 
along the rows and times are arranged along the columns, then a Nx T transmitted data 
matrix X can be written as 



X = 



X(U) - X{\J) 



X{N,\) X{NJ) 



(1) 



Y, a Nx T matrix of received symbols, and V, a A^x T noise matrix, are defined in an 
analogous manner. 

In order to implement one embodiment of the present invention, the present 
inventor designed aNxN equalizer, G, such that 



G(n, :)= argmin4x(«,0-G(/2,:)Y(:,0| 2 J 



(2) 



G(n-) 



If the channel matrix H = diag (H(l) . . . #(N)), and if the received data at time t can be 
written as 



Y(:,0 = HX(:,r) + V(:,0, 



(3) 



where the elements of V(:, t) are uncorrected, then the zero forcing equalizer G can be 
recovered as 



G = diag 



r 1 



(4) 



The minimum mean-square error solution also has a diagonal structure: 
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G = diag 



H(l)* 



\ H ( l )\ 



2 9 
+ - 



Pv(l) 



H(Ny 



\H(N)\ 



2 + P,(A0 



P V (N)J 



(5) 



where Px(k) is the power spectral density (PSD) of the input, and Pv(k) is the PSD of 
the noise on the k th subchannel. These arguments are seen to be the standard one-tap 
FEQ's such as described by T. Starr, J. Cioffi and P. Silverman, Understanding Digital 
Subscriber Line Technology, Upper Saddle River, NJ: Prentice Hall, 1999, referenced 
herein before. 

The present inventor recognized, however, when there is correlation in the noise 
matrix V (e.g., from a correlated noise source such as RFI), it seems reasonable to 
consider equalizers G that have nonzero elements in addition to the main diagonal. One 
embodiment can be exemplified by considering the n th subchannel and letting k(n) be the 
indices of the desired nonzero columns of the n th row of G. Most preferably, G(«, k(n)) 
can then be found such that 



X(n,f) = Y(k(„), t) T G(n, k(n)f 



(6) 



Letting g(w) = G(«, k(n)) r , collecting data from time t = 1, . . ., T, and making the obvious 
substitutions then yields: 



~X{n,\Y 




"Y(k(n),l) r " 


_X(n,T) 




_Y(k(«),r) r 



g(»), 



(7) 



where 



x(n) = Y(«)g(/i) 
Therefore, the n th row of the equalizer G can be recovered as 



(8) 
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g(n) = Y(»rx(it) 



(9) 



This procedure is then repeated for the remaining subchannels to form the complete 
multichannel FEQ: 



G = 



g(l) r 
g(N) T 



(10) 



The present invention is not so limited however, and it shall be understood that 
other solutions, including but not limited to, iterative solutions are also possible using 
least mean-squares (LMS), recursive least-squares (RLS), and their related variants such 
as discussed by G.-O. Glentis, K. Berberidis and S. Theodoridis, "Efficient least squares 
adaptive algorithms for FIR transversal filtering," IEEE Signal Processing Magazine, 
vol. 16, no. 4, pp. 13-41, 1999. The n' h row, for example, of the multichannel FEQ can 
be updated as: 



e(t) = X(n,t)-Y(k(n),t) T g( n ) 



(11) 



g(") = g(>0 + M'M0Y(k(«),0* 



(12) 



for t = 1, . . ., T, where n(t) controls the adaptation. Popular choices include: 



Standard LMS: 



MO = a , 



(13) 



where 
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0<a< — - — (14) 

X, max(R) 



and 



R = £[Y(k(«),0Y(k(«),0"J (15) 
Normalized LMS (NLMS): 



M{t) y9 + Y(k(«),0 // Y(k(«),0' (16) 



where a e (0,2) and 0 < . 

• Power Normalized LMS (PNLMS): 



2 2 

where a 2 (t) = ca 2 (t -Y) + \e(t)\ , cs (0,1), and 0<a< — 

M 



th 

Most preferably, when choosing the nonzero columns of the n row of G, other 
than the n th column (which should be selected as this is the column directly correlated 
with the input through the channel), columns corresponding to subchannels with strong 
correlations with the n th subchannel should be chosen. These are often neighboring 
subchannels, or subchannels corresponding to strong RFI sources. 

Subchannels with strong RFI can be found via any well-known peak picking 
algorithm applied to the noise vector V(:, t). It can be appreciated that when RFI 
dominates the received signal on a subchannel, then adding a weighed version of the 
received subchannel with RFI to the desired subchannel n is (almost) equivalent to the 
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same operation when no data was transmitted on that subchannel. As the RFI and signal 
power become closer, however, not transmitting data on a subchannel allows that 
subchannel to be used more effectively with other subchannels in the multichannel FEQ. 

A good model for the RFI (or other source of subchannel noise correlation) could 
potentially be exploited to improve the performance such as discussed by B. Wiese and J. 
Bingham, "Digital radio frequency cancellation for DMT VDSL," ANSI Tl EL 4/97-460, 
Dec. 1997. A positive aspect of the present method however, is that its adaptive nature 
allows it to be somewhat decoupled from the exact RFI model. Further, the single 
channel FEQ or an existing RFI model could be used to initialize the present adaptive 
process. 

Table 1 below highlights the channel capacity in Mbps for three different 
scenarios including AWGN with no RFI using a single channel FEQ, AWGN with RFI 
using a single channel FEQ, and AWGN with RFI using the present multichannel FEQ. 
The results were obtained by simulating a downstream ADSL channel with RFI added at 
the receiver. Three subchannels were used for the design of the multichannel FEQ for the 
n th subchannel: the n th subchannel, and the two other subchannels with the strongest noise 
power. The multichannel FEQ coefficients were solved for using the NLMS algorithm. 
The channel was modeled as 6000 feet of 26 AWG twisted pair copper wire with -140 
dBm/Hz additive white Gaussian noise (AWGN) and RFI model 1 impairments such as 
discussed by ITU Telecommunications Standardization Sector, G.testbis. Study Group 4, 
Question 15, 2000, referenced herein before. The capacity (Q in bits/s (bps) was 
calculated from: 

C = /olog2(l + £ATC), 
where/c is 4000, and SNR is the gap adjusted signal-to-noise ratio after the FEQ. 
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Table 1 : Capacity calculations 



Noise 


FEQ 


Capacity (Mbps) 


AWGN 


Single channel 


10.951 


AWGN+RFI 


Single channel 


5.794 


AWGN+RFI 


Multichannel 


7.681 



Figure 1 is a graph illustrating a simulated subchannel SNR for 1) additive white 
Gaussian noise (AWGN) with no RFI using a single channel FEQ, 2) AWGN with RFI 
using a single channel FEQ, and 3) AWGN with RFI using the present multichannel FEQ 
in accordance with one embodiment of the present invention. The results suggest that 
significant SNR improvements (and corresponding capacity gains) are possible with the 
present multichannel FEQ. 

Figure 2 is a simplified block diagram depicting an ADSL communication system 
10 using a multichannel FEQ 12 according to one embodiment of the present invention. 
At the transmitter, input data 14 to be transmitted is segmented (serial-to-parallel (S/P) 
converted), as depicted in element 16, into blocks of length N. The n th element is referred 
to as the n th subchannel. The data is transformed (e.g., in the case of ADSL, by an IFFT 
as depicted in element 18, and the addition of a cyclic prefix as depicted in element 20) to 
length iV + P. The transformed data is next unsegmented, as depicted by element 22 
(parallel-to-serial (P/S) converted). The data is then sent through the channel as shown 
by element 24 (perhaps after additional filtering, mixing, etc. not shown). This channel 
data is then received by the ADSL receiver, perhaps after filtering, mixing, etc., and 
generally includes time domain equalization (TEQ) shown by element 26, that is 
conventionally implemented as an FIR filter. Noise 28, such as RFI discussed herein 
before, may be introduced into the signal before it is received by the ADSL receiver. The 
noise that results from RFI at the receiver is deterministically spread by the FFT to 
neighboring subchannels. As a result, instead of one (or a few) subchannels being 
effected by the narrowband impairment, many subchannels are effected. The resulting 
loss of signal-to-noise ratio (SNR) can unacceptably limit the capacity of the 
communication channel 24. With continued reference now to Figure 2, the data is S/P 
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converted as shown by element 30, into blocks of length Af+ P, subsequent to time 
domain equalization 26. The data is then transformed (e.g., in the case of ADSL, by 
removing the cyclic prefix as shown by element 32 and a FFT as shown by element 34) 
to length N. Following this transformation process, the present multichannel FEQ 12 is 

th 

applied to each block of data. The multichannel FEQ 12 output data for the n 
subchannel can then be expressed as 



x(nJ) = Y(k(n)j) T g(n). 



as discussed herein before with reference to equation (7). The output data from the 
multichannel FEQ 12 is the P/S converted, as shown by element 36, to form an estimate 
of the original data. 

Figure 3 is a flow diagram illustrating a method 100 of estimating multichannel 
FEQ coefficients for the n th subchannel, g(w), according to one embodiment of the present 
invention. As discussed above, the present method relates to T blocks of transmitted 
data, X(:, t),t = 1, . . T, and Tblocks of received data, Y(:, t), t = 1, . . r depicted by 
elements 102 and 104 respectively. The received data can then be used to determine a set 
of subchannels, k(«), for the multichannel FEQ 12 to use in conjunction with subchannel 
h, for each subchannel n used to transmit data, as shown by element 106. The set of 
subchannels, k(w), includes the n th subchannel and one or more of the following: the n th 
subchannel's neighboring subchannels; subchannels where RFI is located; and 
subchannels where the noise power is (relatively) strong. Finally, the multichannel FEQ 
coefficients for the n th subchannel, g(«), are generated as shown by element 108, via 
equation (2) using all of the received data at once, according to equation (9), or 
adaptively, using a block of data at a time, according to equations (1 1) and (12). These 
multichannel FEQ coefficients are then used in a complex multiply operation to adjust for 
the gain and phase of the effective channel response in each subsymbol. 

In view of the above, it can be seen the present invention presents a significant 
advancement in the art of communication between subscribers over shared media. 
Further, this invention has been described in considerable detail in order to provide those 
skilled in the data communication art with the information needed to apply the novel 
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principles and to construct and use such specialized components as are required. In view 
of the foregoing descriptions, it should be apparent that the present invention represents a 
significant departure from the prior art in construction and operation. However, while 
particular embodiments of the present invention have been described herein in detail, it is 
to be understood that various alterations, modifications and substitutions can be made 
therein without departing in any way from the spirit and scope of the present invention, as 
defined in the claims which follow. 
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